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ABSTRACT 

The protostellar jets driven by the formation of the first stars are studied 
by using three-dimensional magnetohydrodynamical (MHD) nested grid simula- 
tions. Starting from a slowly rotating spherical cloud of 5.1 x 1O 4 M permeated 
by a uniform magnetic field, we follow the evolution from the central number den- 
sity n c = 10 3 cm~ 3 (where the radius of the object r = 6.6 pc) to n c ~ 10 23 cm -3 
( r ~ 1-R©)- We resolve the cloud structure more than 8 orders of magnitude 
in spatial extent and 20 orders in density contrast. We calculate four mod- 
els that differ in initial magnetic field strengths and angular velocities. In all 
models, protostars of ~ 10~ 3 M & are formed at n c ~ 10 22 cm -3 in accordance 
with one-dimensional calculations. By this epoch, the magnetic flux density is 
amplified by 10 orders of magnitude from the initial value. Consequently, the 
formed protostar possesses the magnetic field of ~ 10 6 G that is much larger 
than the flux density of the present counterparts, reflecting the fact that the 
dissipation of a magnetic field is ineffective in primordial gas clouds. If the initial 
magnetic field B > 10~ 9 (n c /10 3 cm~ 3 ) 2 / 3 G, the protostellar jet is launched and 
its velocities reaches ~ 70 km s _1 by the time the protostellar mass becomes 
(4 — 6) x 1O~ 3 M , and a fraction (3 — 10%) of the accreting matter is blown 
off from the central region. Owing to the interaction of these ejecta with sur- 
rounding matter, expanding bow shocks are created at both heads of the jet. If 
this jet continues to sweep out the surrounding gas that otherwise accretes onto 
the central star or circumstellar disk, the final mass of the first star can be sub- 
stantially reduced. In addition, dense post-shock regions behind the bow shocks 
are expected to promote the chemical reactions (formation of H 2 and HD), and 
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this provides possible environments for subsequent low-mass star formation in 
the early universe. 

Subject headings: cosmology: theory — early universe — galaxies: formation — 
ISM: jets and outflows — MHD — stars: formation 

1. Introduction 

Magnetic field, which is an indispensable ingredient in the present-day local interstellar 
medium (ISM), has been largely overlooked in studies of the first star formation. This is 
probably because the magnetic force has been believed to be unimportant in comparison 
with either the gravity or the thermal pressure in the early universe. 

Theoretical studies showed that a collapsing primordial gas follows a self-similar solu- 
tion with pressure P in the cooling gas scaling with density p as P oc p 7cff where 7 c fr ~ 1.1 
(e.g., Omukai & Nishi 1998). In a gravitationally collapsing state such as in self-similar 
solutions, the relative importance of pressure and gravity inside of the runaway collapsing 
core remains almost the same. In contrast, the magnetic field in the homologously collapsing 
medium scales as B ~ p 2 / 3 , and thus, the magnetic pressure increases more rapidly than the 
gas pressure. In general, if the magnetic field scales as B ~ p 7s and 7# > 7 c fr 5 the dynamical 
importance of the magnetic field increases continuously in the collapsing primordial cloud. 
Therefore we cannot entirely neglect the effect of the magnetic field in the collapsing primor- 
dial gas even if it is relatively small initially. Ichiki et al. (2006) showed that magnetic fields 
are created by the cosmological fluctuations when they enter the cosmic horizon before the 
cosmic recombination epoch. By this mechanism, the magnetic field of ~ 10~ 15 — 10~ 14 G is 
produced at the scale of the first cosmological objects (~ 10 5 ~ 6 M Q ). Langer et al. (2003) 
proposed another mechanism for magnetic field generation during the epoch of the cosmic 
reionization. They found that magnetic fields in the intergalactic matter on the order of 
~ 10~ n G are created at the scale of protogalaxies, which increase up to ~ 10~ 7 — 10~ 8 G 
if the objects experience the spherical collapse up to number density n = 10 3 cm~ 3 . These 
fields can influence the evolution of primordial gas. 

In this paper, we study the evolution of weakly magnetized primordial star-forming cores 
and the formation of protostars using three-dimensional numerical magnetohydrodynamical 
(MHD) simulations. As a result, we demonstrate that protostellar jets are driven by the 
formed protostars even in the primordial condition. 
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2. Model and Numerical Method 

To study the evolution of primordial star-forming cores taking account of a magnetic 
field and self-gravity in a large dynamic range of density and spatial scale, we use the 
three-dimensional ideal MHD nested grid method. The ideal MHD equations including the 
self-gravity are 

f + V-( P «) = 0, (1) 

dv I 
p— + p(v • V)« = -VP - — B x (V x B) - pV0, (2) 
ot 4tt 

°— = V x (v x B), (3) 
V 2 = AnGp, (4) 

where p, v, P, B, and denote the density, velocity, pressure, magnetic flux density, 
and gravitational potential, respectively. The ideal MHD approximation is justified in 
our case. Studying the coupling of the magnetic field in star-forming cores, Maki & Susa 
(2004) concluded that magnetic field is frozen to the gas for the initial field strength below 
10~ 5 (nH/10 3 cm~ 3 )°- 55 G. For gas pressure, we use a barotropic relation that approximates 
the result of one-zone calculation where thermal and chemical processes for primordial gas 
are solved in detail (Omukai et al. 2005). This relation (thick solid line) as well as the 
original data (thick dotted line) for pressure evolution is plotted as a function of number 
density in the upper panel of Figure 1. To stress variations of pressure with density, we plot 
P/n, which is proportional to the gas temperature if the mean molecular weight is constant. 

As the initial states, we take a spherical cloud whose density is twice higher than that for 
hydrostatic equilibrium with external pressure (i.e., so called Bonnor-Ebert sphere, Bonnor 
1956; Ebert 1955). The cloud rotates rigidly (fi ) around the ^-axis and is threaded by 
a uniform magnetic field Bq parallel to the rotation axis. The initial central density is 
taken as n c = 10 3 cm -3 , which corresponds to the epoch of formation of star-forming cores 
(Abel et al. 2000; Bromm et al. 1999, 2002). The radius and mass inside the sphere 
are R c = 6.6 pc and M c = 5.1 x 1O 4 M , respectively. The initial temperature is 250 K. 
After the formation of the star-forming core, self-gravity of gas dominates the dark matter 
gravity (Abel et al. 2002). Thus, we do not include the dark matter component in our 
calculation. The models are characterized by two parameters: initial magnetic field strength 
(B ) and initial angular velocity (fio). We calculate four models denoted as B6, B7, B8 and 
B9, whose model parameters are (B , Q ) = (10~ 6 G, 3.3 x 10~ 16 s _1 ; Model B6), (10~ 7 G, 
3.3 x 10~ 17 s _1 ; Model B7), (10~ 8 G, 2.1 x 10~ 17 s _1 ; Model B8), and (10~ 9 G, 1.3 x 10~ 17 
s _1 ; Model B9), respectively. We have observed fragmentation of clouds in cases with faster 
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rotation. For example, with initial magnetic field B = 1CT 6 G, the cloud fragments if 
angular velocity Q > 3.3 x 1CT 16 s -1 . Our rotation rates are close to the critical values 
for fragmentation according to our numerical experiment. Since our focus in this paper is 
on driving of protostellar jet rather than on core fragmentation, we adopt somewhat slow 
rotation rates in comparison with some examples of Bromm et al. (2002) who studied the 
fragmentation of primordial clouds by hydrodynamical simulation. We leave the stage of the 
evolution with larger rotation rate to our subsequent papers. 

We adopt the nested grid method (for detail, see Machida et al. 2005a,b) to achieve 
high spatial resolution near the origin. Each level of a rectangular grid has the same number 
of cells (= 64 x 64 x 32), although the cell width h(l) depends on the grid level /. The 
cell width is reduced by a factor of two for every upper level. The calculation is started 
with three grid levels (/ = 1,2,3). Box size of the coarsest grid / = 1 is chosen to be 
2R C . At the boundary, the magnetic field and ambient gas are assumed to rotate with the 
constant angular velocity f2 (for detail see Matsumoto & Tomisaka 2004). A new finer grid 
is generated whenever the minimum local Jeans length Aj becomes smaller than 8/i(/ max ), 
where h is the cell width. The maximum level of grids is restricted to Z max ^ 30. Since the 
density is highest in the finest grid, generation of new grid ensures the Jeans condition of 
Truelove et al. (1997) with a margin of a safety factor of 2. 



3. Results 

Starting from the number density n c = 10 3 cm~ 3 , the evolution of primordial clouds 
is calculated up to ~ 10 23 cm~ 3 . Figure 1 shows the evolution of central angular velocities 
(f2 c ; upper panel) and magnetic flux densities (B c ; lower panel) against the central number 
density (n c ) for all models. Both of angular velocity and magnetic flux density in model 
B6 increase in proportion to nj 3 , which is the case of spherical collapse, for 10 3 cm~ 3 < 
n c ^ 10 14 cm~ 3 (see Machida et al. 2005a, 2006a). For higher densities, the growth rates 
in unit of logarithmic density (d\ogB c /d\ogn c and d\ogQ c /d\ogn c ) become smaller, and 

1/2 

the angular velocity and magnetic flux density increase in proportion to n c , because cloud 
collapses non-homologously into the self-similar disk-like structure (Machida et al. 2005a, 
2006a). On the other hand, in models B7, B8, and B9, both the magnetic flux densities 
and angular velocities continue to increase in proportion to n^ 3 until n c ~ 10 21 cm~ 3 . This 
indicates that, owing to small magnetic field and slow rotation, disk formation is delayed and 
spherical collapse continues until very high density. In all models, the angular velocities have 
converged to the same scaling law before n c ~ 10 19 cm~ 3 owing to the change of geometry 
of the collapse and magnetic braking: for stronger magnetic field, the magnetic braking is 
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more effective and the growth rate of angular velocity tends to be smaller. Since we assume 
rapid (slowly) rotation for strongly (weakly, respectively) magnetized clouds in this paper, 
the angular velocities tend to converge to the unique scaling law as the collapse proceeds. 

At n c ~ 10 20 cm~ 3 that corresponds to the completion of molecular hydrogen dissoci- 
ation, the equation of state becomes adiabatic as shown in the upper panel of Figure 1. 
The hydrostatic core of mass M* ~ 10~ 3 M Q is formed at n c ~ 10 22 cm~ 3 in all models in 
agreement with one-dimensional calculation (Omukai & Nishi 1998). We call this hydro- 
static core a protostar hereafter. The protostar rotates at angular velocity of Q ~ 10~ 4 s -1 
(~ 8.63 day -1 ) initially, which is comparable to or slightly faster than those of the present- 
day pre-main sequence stars (Herbst et al. 2002). By this time, the magnetic field has been 
amplified to the level of B c ~ 6 x 10 5 G (Model B9) - 6.1 x 10 6 G (Model B6) that is much 
stronger than that of the present day protostar or, so called, the second core (see Machida 
et al. 2006b). This is because the magnetic field in primordial gas does not dissipate either 
by Ohmic dissipation or ambipolar diffusion, as shown by Maki & Susa (2004). The formed 
protostars in the present-day star formation process have the magnetic field of only ~ 10 3 
G, owing to effective Ohmic dissipation (Machida et al. 2006b). 

Figure 2 shows the density distributions and velocity vectors on the y = cut planes for 
models B7, B8, and B9 at the end of the calculation. Each upper panel is a close-up view of 
the corresponding lower panel. The level of the outermost grid is indicated at the upper left 
corner of each panel. Finer grids are overplotted by boxes in each panel. We omit model B6 
from this figure, since the evolution of model B6 is almost the same as that of model B7. 
We can see in Figure 2a that a strong jet is driven from the central protostar. 

A spherical protostar of radius 0.65 R Q forms at the center, which is surrounded by 
a disk-like structure extending up to ~ 2R Q . The jet starts flowing 49.1 hours after the 
protostar formation epoch. 1 Rippling of density contours in Figure 2 indicate the occurrence 
of strong mass ejection near the center. The mass of outflowing gas is 7 x 10~ 4 M , while the 
mass of the protostar and the disk are 4.4 x 10~ 3 M and 1.6 x 10~ 3 M , respectively, at the 
end of the calculation. Namely, about 10% of the accreted matter is ejected as the jet. The 
bow shocks are generated at the top of the jet, z ~ ±0.1 AU (Fig. 26), owing to interaction 
with ambient matter. The jet is well-collimated and creates a cocoon-like structure as shown 
in Figure 2b. At this epoch, the maximum speed of the jet is 66.3 km s -1 . In model B6, the 
mass accretion rate, outflow rate, the speed of the jet, and density structure are very similar 
to those in model B7. 



1 We define the protostar formation epoch as the time when the central density reaches n c = 10 20 cm 
where the equation of state becomes adiabatic. 
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In model B8, the jet appears 45 hours after the protostar formation epoch. The speed 
of the jet in model B8 is slower than that in model B7. We can see from the arrows in 
Figure 2 c and d that the speeds of jet are comparable to or slower than the infall speed in 
model B8, while the speeds of the jet in model B7 are much larger than the infall speed 
(Fig. 2a and b). However, the rippling density contours indicative of strong mass ejection 
are also seen in this model (Fig. 2 c). The mass of the core, disk, and outflowing gas increase 
up to 6.15 x l(r 3 M , 1.25 x 1(T 3 M , and 2.5 x 1O~ 4 M , respectively, by the end of the 
calculation. Thus, about 3% of the accreted matter is converted to outflow until the end 
of the calculation in this model. Owing to smaller magnetic field, the outflow efficiency is 
reduced by about a factor of three from model B7. 

Further reduction of magnetic field results in the absence of jet in model B9. In this 
model, an oblate protostar forms inside a thin disk as shown in Figure 2e and /. Although 
we followed the evolution of the accreting protostar up to 201 hours after the protostar 
formation epoch, we could not observe any sign of the jet. Thus, without hindrance by jet, 
the gas can accrete onto the protostar directly in the vertical direction as well as through the 
disk, while, in other models (B6-8), the gas accretes mostly through the thin disk. After the 
central protostar forms, magnetic field inside the star is significantly amplified both by the 
magneto-rotational instability (Balbus & Hawley 1991) and by the shearing motion between 
the star and the ambient medium. Consequently, the plasma beta j3 p , which is the ratio of 
gas pressure to magnetic pressure, inside the star becomes almost unity, i.e. the magnetic 
energy is comparable to the thermal energy. Just outside the disk, however, the magnetic 
energy is much smaller than the thermal energy (/3 p 1). This small magnetic energy 
compared with the thermal energy is the reason why the jet is not launched in model B9. 
In contrast, in models B6-8, the magnetic energy is comparable to or lager than the thermal 
energy even outside the disk. Note, however, that a long time calculation of model B9 remain 
to be carried out to clarify whether a jet is eventually launched or not. 



4. Conclusion and Discussion 

From these experiments, we conclude that the protostellar jet is driven by the primordial 
protostar if the initial magnetic field B > 10~ 9 G at n — 10 3 cm~ 3 . This critical value 
corresponds to the magnetic field of 5 x 10~ 13 G in the background material (n c ~ 0.01 cm -3 ) 
if the spherical collapse (B oc p 2 ^ 3 ) is assumed. Whether this initial magnetic field is realized 
or not remains uncertain. The magnetic field expected from the cosmological fluctuations 
falls below this value about an order of magnitude (Takahashi et al. 2005; Ichiki et al. 
2006). However, with slightly more amplification by an unknown mechanism, the field could 
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exceed the critical value. The field created at the cosmic reionization (~ 10~ n G) seems to 
be sufficient (Langer et al. 2003). However, after the reionization epoch, the primordial 
star formation occurs only in large halos with the virial temperature exceeding 10 4 K, being 
strongly suppressed in smaller halos by UV background radiation (e.g., Kitayama et al. 
2001). Even if the initial magnetic field is weaker than our critical value, the protostellar jet 
could be launched in a later phase of star formation. By analytically studying the evolution 
of accretion disks around the first stars, Tan & Blackman (2004) suggested that magnetic 
fields amplified in the disk eventually give rise to protostellar outflows during the protostellar 
accretion phase. Since we did not follow the evolution up to the late phase considered by 
Tan & Blackman (2004), we cannot draw any conclusions concerning their expectation. 

Note that we assumed rather slow rotation of clouds (f2 = 3.3xl0 _16 s _1 — 1.3xl0 _17 s _1 ) 
as the initial state. Strictly speaking, this critical value for the magnetic field strength is 
only valid for clouds having small angular velocities. When an initial cloud rotates rapidly, 
the jet might not be driven by the protostar because the amplification rate of magnetic field 
changes according to the geometry of the collapse and cloud fragmentation (Machida et al. 
2005a, 2006a). This possibility will be studied in our subsequent works. 

Although numerical studies on fragmentation of first cosmological objects concluded 
that fragmentation mass scale of metal-free gas is as massive as 100 — 1000 M© (e.g., Bromm 
et al. 1999, 2002; Abel et al. 2000; Omukai & Palla 2003), they are limited to the case of 
star formation from the cosmological initial condition, which is unaffected by astrophysical 
phenomena. Jet and outflow from the first stars disturb the condition of star forming regions 
in the early universe. What kinds of stars are formed in such regions? Observations of local 
star-forming regions indicate that the protostellar outflows compress surrounding matter and 
induce subsequent star formation (e.g., Bally et al. 1999; Quillen et al. 2005). Since our 
calculations demonstrated that the jets are indeed driven by the first stars, similar events 
may possibly occur also in the case of the first stars. 

When the velocity of jets exceeds ~ 50 km s _1 , gas is heated above ~ 10 5 K and then 
ionized. If the shocked layer becomes fully ionized, a large amount of molecular hydro- 
gen forms in the postshock flow by using electron as a catalyst during slow recombination 
(MacLow & Shull 1986; Shapiro k Kang 1987; Susa et al. 1998). The enhanced H 2 cooling 
realizes low-temperature environment (< 150K) favorable for HD formation, and this HD 
cooling further reduces the temperature below ~ 100 K. In such cold clumps, low-mass stars 
are expected to be produced (e.g., Uehara & Inutsuka 2000; Nakamura & Umemura 2002; 
Machida et al. 2005c; Nagakura & Omukai 2005). Those low-mass stars survive to the 
present and could be observed in the halo. Suda et al. (2004) indicated that the extremely 
metal-poor ([Fe/H]<-5) stars (Christlieb et al. 2001; Frebel et al. 2005) discovered re- 
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cently were formed as binary members from metal-free gas, and then have been polluted by 
companion stars during the stellar evolution. 

However, there is a caveat in discussing the subsequent star formation around the first 
stars. UV photons emitted by the formed massive star may dissociate H 2 and HD molecules 
in the surrounding matter thereby suppressing the subsequent star formation (Omukai & 
Nishi 1999). Future study on induced star formation by protostellar outflow should include 
this effect. 

Not only the subsequent star formation but also the first star formation process itself 
is affected by the outflow. For example, in the present-day star formation scenario, the final 
stellar mass is envisaged to be affected by the mass outflow (for low-mass stars) or radiative 
feedback (for massive stars) from the protostar. These mechanisms stop mass accretion from 
the parent cloud to the protostar, as well as blow off the parent cloud itself (Nakano et al. 
1995). In the case of primordial star formation, the radiative feedback is insufficient for 
halting the accretion of metal-free gas because of the lack of dust grain (Omukai & Palla 
2001; 2003). The mass of first stars grows up to ~ 600 M Q if the mass accretion continues 
unimpeded. The outflow from the first star may reduce its final stellar mass by blowing off 
a considerable fraction of the cloud matter before accreting onto the star. At the same time, 
the jet may reduce the star formation efficiency by diverting material from the accretion 
disk. 

The jets from the first stars also contribute to disperse the magnetic field lines into the 
intergalactic medium. This may affect the formation of next generation stars and galaxies. 
To study this effect further, we need to calculate the evolution up to the later phase possibly 
by adopting implicit numerical integration method for MHD equations. 

We greatly benefited from discussion with T. Kato, K. Saigo, and H. Susa. We also thank 
T. Hanawa for contribution in developing the nested grid code. Numerical calculations were 
carried out with a Fujitsu VPP5000 at the Astronomical Data Analysis Center, the National 
Astronomical Observatory of Japan. This work was supported partially by the Grants-in-Aid 
from MEXT (15740118, 16077202, 16740115). 



REFERENCES 

Abel, T., Bryan, G. L., & Norman, M. L. 2000, ApJ, 540, 39 
Abel, T., Bryan, G. L., & Norman, M. L. 2002, Science, 295, 93 



- 9- 

Bally, J., Reipurth, B., Lada, C. J., & Billawala, Y. 1999, AJ, 117, 410 
Balbus, S. A., & Hawley, J. F. 1991, ApJ, 376, 214 
Bonnor, W. B. 1956, MNRAS, 116, 351 

Bromm, V., Coppi, P. S., & Larson, R. B. 1999, ApJ, 527, L5 

-, 2002, ApJ, 564, 23 
Christlieb, N., Green, P. J., Wisotzki, L., & Reimers, D. 2001, A&A, 366, 898 
Ebert, R. 1955, Z. Astrophys., 37, 222 
Frebel, A., et al. 2005, Nature, 434, 871 

Herbst, W., Bailer- Jones, C. A. L., Mundt, R., Meisenheimer, K., & Wackermann, R. 2002, 
A&A, 396, 513 

Ichiki, K., Takahashi, K., Ohno, H., Hanayama, H., & Sugiyama, N., 2006, Science, 311, 787 

Kitayama, T., Susa, H., Umemura, M., & Ikeuchi, S. 2001, MNRAS, 326, 1353 

Langer, M., Puget, J.-L., & Aghanim, N. 2003, Phys. Rev. D, 67, 043505 

Machida, M. N., Matsumoto, T., Tomisaka, K., & Hanawa, T. 2005a, MNRAS, 362, 369 

Machida, M. N., Matsumoto, T., Hanawa, T., & Tomisaka, K. 2005b, MNRAS, 362, 382 

Machida, M. N., Tomisaka, K., Nakamura, F., & Fujimoto, M. 2005c, ApJ, 622, 39 

Machida, M. N., Matsumoto, T., Hanawa, T., & Tomisaka, K. 2006a, ApJ accepted 

Machida, M. N., Inutsuka, S., & Matsumoto, T., 2006b, ApJ submitted 

MacLow, M. M., & Shull, J. M., 1986, ApJ, 302, 585 

Maki, H., & Susa, H., 2004, ApJ, 609, 473 

Matsumoto, T., & Tomisaka, K. 2004, ApJ, 616, 266 

Nagakura, T., & Omukai, K. 2005, MNRAS, 364, 1378 

Nakamura, F., & Umemura, M. 2002, ApJ, 569, 549 

Nakano, T., Hasegawa, T., & Norman, C. 1995, ApJ, 450, 183 



Omukai, K. & Nishi, R. 1998, ApJ, 508, 141 
Omukai, K. & Nishi, R. 1999, ApJ, 508, 141 
Omukai, K., & Palla, F. 2001, ApJ, 561, L55 
Omukai, K., & Palla, F. 2003, ApJ, 589, 677 

Omukai, K., Tsuribe, T., Schneider, R., & Ferrara, A. 2005, ApJ, 626, 627 

Quillen, A., Frank, A., Cunningham, A., Blackman, E., & Pipher, J. 2005, ApJ, 631, 941 

Shapiro, P. R., & Kang, H., 1987, ApJ, 318, 32 

Suda, T., Aikawa, M., Machida, M. N., Fujimoto, M. Y., & Iben, I. J. 2004, ApJ, 611, 476 
Susa, H., Uehara, H., Nishi, R., & Yamada, M. 1998, Prog. Theor. Phys., 100, 63 
Takahashi, K., Ichiki, K., Ohno, H., & Hanayama, H. 2005, Phys. Rev. Lett., 95, 121301 
Tan, J. C., & Blackman, E. G. 2004, ApJ, 603, 401 

Truelove J, K., Klein R. L, McKee C. F., Holliman J. H., Howell L. H., & Greenough J. A., 
1997, ApJ, 489, L179 

Uehara, H., & Inutsuka, S. 2000, ApJ, 531, L91 



This preprint was prepared with the AAS IATgX macros v5.2. 



- 11 - 




n (cm" 3 ) 



Fig. 1. — Upper panel: The gas pressure normalized by number density derived from a one- 
dimensional calculation (thick solid line) and our fit used in the numerical simulation (thick 
dotted line) are plotted. The central angular velocities (f2 c ) against the number density (n c ) 
are also plotted for models B6 (solid line), B7 (broken line), B8 (dotted line), and B9 (dash- 
dotted line). Lower panel: The magnetic flux densities (B c ) against the central density for 
models B6 (solid line), B7 (broken line), B8 (dotted line), and B9 (dash-dotted line). The 
relation oc p 2 / 3 , valid for the spherical collapse, is also plotted for comparison. 



Fig. 2. — The density (color-scale and white contours) and velocity distribution (arrows) 
on the cross-section in the y = plane for models B7 (panels [a] and [b]), B8 (panels [c] 
and [d]), and B9 (panels [e] and [f]). Model name, elapsed time, central number density, 
and velocity scale are plotted in each upper panel. Grid level (/) are displayed at upper left 
corner. Each upper panel is a close-up view of the corresponding lower panel. 



